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ABSTRACT 

We present new quasars discovered in the vicinity of the Andromeda and 
Triangulum galaxies with the LAMOST (Large Sky Area Multi-Object Fiber 
Spectroscopic Telescope, also named Guoshoujing Telescope) during the 2010 and 
2011 observational seasons. Quasar candidates are selected based on the available 
Sloan Digital Sky Survey (SDSS), Kitt Peak National Observatory (KPNO) 4 m 
telescope, Xuyi Schmidt Telescope Photometric Survey (XSTPS) optical, and 
Wide-field Infrared Survey Explorer (WISE) near infrared photometric data. We 
present 509 new quasars discovered in a stripe of ~135 deg^ from M 31 to M 33 
along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new 
quasars discovered in an area of ~ 100 deg^ that covers the central region and the 
southeastern halo of M31 in the 2010 commissioning datasets. These 526 new 
quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. 
They represent a significant increase of the number of identified quasars in the 
vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this 
region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, 
of which 5, 20 and 75 are newly-discovered. These bright quasars provide an 
invaluable collection with which to probe the kinematics and chemistry of the 
interstellar/intergalactic medium (ISM/IGM) in the Local Group of galaxies. A 
total of 93 quasars are now known with locations within 2.5° of M 31, of which 73 
are newly discovered. Tens of quasars are now known to be located behind the 
Giant Stellar Stream, and hundreds behind the extended halo and its associated 
substructures of M 31. The much enlarged sample of known quasars in the vicinity 
of M31 and M33 can potentially be utilized to construct a perfect astrometric 
reference frame to measure the minute proper motions (PMs) of M31 and M33, 
along with the PMs of substructures associated with the Local Group of galaxies. 
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Those PMs are some of the most fundamental properties of the Local Group. 

Subject headings: galaxies: individual (M31, M33) — quasars: general — quasars: 
emission lines 
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1. Introduction 

Being the most luminous member of the Local Group and the nearest archetypical 
spiral galaxy, M31 serves as the best astrophysical laboratory for the studies of the physical 
and astrophysical processes that govern the structure, kinematics and the formation 
and evolution of distant galaxies. Recent deep optical surveys have revealed complex 
substructures within hundreds of kiloparsec (kpc) of M31, with some of them stretching 
from M31 all the way to M33, which is about 200 kpc from M31, suggesting a possible 
recent close encounter of the two galaxies (Ibata et al. 2007; McConnachie et al. 2009). 
Detailed chemical and kinematic investigations of M31 and associated substructures are 
vital for our understanding of M31, and also for the theory of galaxy formation and 
evolution in general. 

Finding background quasars in the vicinity of M31 have the following two potential 
important applications. Firstly, the background quasars can serve as an ideal reference 
frame which would allow highly accurate astrometric measurement of the minute proper 
motion (PM) of M31 and also the PMs of its associated coherent substructures. Secondly, 
absorption-line spectroscopy of bright background quasars can be used to probe the 
distribution, chemical composition and kinematics of the interstellar medium (ISM) of M31, 
the Milky Way and the intergalactic medium (IGM) of the Local Group of galaxies. This 
technique can be used to probe structures further out from the center of the galaxy than 
would be possible with traditional stellar spectroscopy, or using H I 21cm observations. In 
addition, M 31 has alway been a focus of interest for astronomical observations. Extensive 
photometric data spanning tens of years exist. Therefore, identifying background quasars 
in the vicinity of M31 should also enable the construction of a valuable sample for studying 
quasar variability. 

The PM of M31 represents one of the most fundamental quantities of the Local Group, 



-6- 

and its accurate measurement is vital to the study of the formation, evolution and dynamics 
of the Local Group, in particular for quantifying the dark matter content. M 31 is known to 
be located at a distance of ~ 785 kpc (McConnachie et al. 2005), and is moving towards the 
Milky Way with a velocity of 117 km s^^ (Binney & Tremaine, 1987, p. 605). However, its 
transverse velocity has remained unknown for a long time. The Very Long Baseline Array, 
with a PM measurement uncertainty of the order of one micro-arcsec, has measured the 
distance and PM of M33 by water maser observations (Brunthaler et al. 2005). However, 
water maser sources in M31 have been discovered only recently (Darling 2011), thus we do 
not yet have a long enough time baseline for a meaningful PM measurements. Based on 
kinematic information about M31's satellite galaxies, M33 and IC 10, Loeb et al. (2005) 
and van der Marcel & Guhathakurta (2008) deduced theoretical estimates for the PM of 
M31 of about 80 km s~^ (20 /las yr^^). Sohn, Anderson & van der Marel (2012) presented 
the first direct PM measurements for three M31 fields (~ 2.7x2.7 sq. arcmin. for each 
field) using HST observations that span a time baseline of 5-7 yrs. Due to the sparseness 
of known background quasars, their analysis had to resort to using compact background 
galaxies as reference sources, even though it was hard to obtain highly accurate positions for 
those galaxies given their nature as extended sources. Being essentially point sources with 
zero PMs, background quasars can serve as perfect reference sources for highly accurate PM 
measurements, provided that enough of them with sufficient space density can be identified. 

At the distance of M31, even a luminous red giant branch star has an / magnitude 
fainter than 20, hence high resolution spectroscopic determinations of their chemical 
composition are no easy tasks even for a 10 m class telescope (Ibata et al. 2005; Gilbert 
et al. 2009). Absorption-line spectroscopy of background quasars allow one to probe the 
distribution, chemical composition and kinematics of the ISM associated with M31, and this 
technique can be used to probe material further out from the galaxy center than is possible 
with traditional stellar spectroscopy or with H I 21cm observations (Chemin, Carignan & 
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Foster 2009). A number of papers have been published (see Savage et al. 2000; Schneider et 
al. 1993 and references therein), reporting Milky Way absorption line systems which were 
detected in quasar spectra. Analyses of those systems have revealed great complexity of the 
Milky Way halo gas, which exhibits a wide range of ionization states, chemical compositions 
and kinematics. Similar absorption- line spectroscopic analysis of M31 background quasars 
will help understand the formation history of M 31 and the Local Group. 

A number of efforts have been taken to search for background quasars behind nearby 
galaxies, with the intended purpose being to study their PMs or ISM/IGM. Tinney, Da 
Costa & Zinnecker (1997) and Tinney (1999) obtained a sample of quasars behind the 
Milky Way satelite galaxies. Kim et al. (2012) published 663 background quasar candidates 
in the Large Magellanic Cloud (LMC), Kozlowski et al. (2012; see also references therein) 
increased the number of known background quasars to 200 behind the southern edge of the 
LMC. For M31, Crampton et al. (1997) presented several quasar candidates. Untill recently, 
Sloan Digital Sky Survey (SDSS; York et al. 2000) obtained 75 low-redshift quasars in three 
spectroscopic plates in two fields in the outer halo of M31 (c.f. Adelman-McCarthy et al. 
2006, 2007). As far as we are aware, since 2007 no further observations targeting quasars 
around M31 have been made by SDSS. In the first of this series of papers we presented 14 
new quasars discovered with the LAMOST during its early commissioning phase (see Huo 
et al. 2010, hereafter Paper I). 

LAMOST is a quasi-meridian reflecting Schmidt telescope with an effective aperture 
of about 4m and a field of view of 5deg in diameter (Wang et al. 1996; Su et al. 1998; 
Xing et al. 1998; Zhao 2000; Cui et al. 2004; Zhu et al. 2006; Cui et al. 2010, and 



Cui, Zhao & Chu et al. 2012, c.f. http://www.lamost.org/website/en/). Recording 4000 



celestial object spectra simultaneously with a parallel fully controllable fiber positioning 
system, LAMOST is the telescope with the highest rate of spectral acquisition at the 



present time. After a period of two years of commissioning for the purpose of technical 
fine-tuning, performance characterization and optimization, as well as scientific capability 
demonstration, the LAMOST pilot surveys were initiated in October 2011 and completed in 
spring the following year (see Zhao et al. 2012). In this paper, we report newly-discovered 
background quasars in the vicinity of M31 and M33, based on LAMOST data obtained 
during the commissioning phase of 2010 as well as during the pilot surveys from autumn 
2011 to spring 2012. 

2. Candidate Selection 

With the aim of searching for background quasars in the vicinity of M31 and M33, we 
selected quasar candidates based on the available optical, near infrared (IR) photometric 
data in the extended area of M31 and M33. Only point-source quasar candidates were 
selected in order to provide a clean sample for the PM measurements and ISM/IGM studies. 

Fig. [1] shows the fields observed with the LAMOST during 2010 and 2011 observational 
seasons. M31 is covered by a special scan of a total area of 100 square degrees, 5° wide and 
20° long, along the northeast-southwest direction (Adelman-McCarthy et al. 2006, 2007). 
M33 is covered by the SDSS southern Galactic cap survey (Aihara et al. 2011). Following 
Richards et al. (2002), low-redshift quasar candidates are selected as outliers from the 
stellar loci on the SDSS color-color diagrams (see also Paper I for the color-color diagrams 
and selection criteria). The candidates are required to have magnitude errors less than 
0.2 magnitude in ugriz. With a magnitude limit set at « < 20.0, the low-redshift quasar 
candidates thus selected have a number density of 27/sq. deg. 

As shown in Fig.[Tl there are no reliable photometric data from the SDSS in the central 
area of M31 and M33 given the very high stellar densities there. To select candidates in 
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those crowded fields, we have made use of the UBVRI photometric catalog of M31 and 
M 33 provided by the Kitt Peak National Observatory (KPNO) 4m telescope (Massey et al. 
2006). In doing so, we convert the KPNO UBVRI magnitudes to the SDSS ugriz using 
the transformations provided by Jester et al. (2005) for quasars of redshifts z < 2.1. 

In addition to those purely optically-selected candidates, we have selected quasar 
candidates combining the optical and IR photometry. Sources from Wide-field Infrared 
Survey Explorer (WISE; Wright et al. 2010) with colors Wl — W2 > 0.8 were selected 
(see Fig. 11, 12 in Yan et al. 2013 for the color-color distributions of different source 
classification; also Stern et al. 2012; Wu et al. 2012 and references therein) and 
cross-correlated with the optical point-source catalogs of SDSS, KPNO and the Xuyi 
Schmidt Telescope Photometric Survey (XSTPS; Liu et al., in preparation, see also 



http://kiaa.pku.edu.cn/DSSGAC/overview.html). XSTPS has a component that extends 



to the M31/M33 area, as shown in Fig.[Tl In total, about 26 quasar candidates per sq. 
degree were selected by this technique for a limiting magnitude i < 20.0. For the sky area 
also covered by the SDSS, we find that about 60% of these WISE IR quasar candidates 
overlap with the optically-selected low-redshift sample. As shown in Wu et al. (2012), this 
IR-based technique is capable of finding quasars with redshifts up to 2 < 3.5. 

On the other hand, for the central regions of M31 and M33, this IR selection algorithm 
also yields too many interlopers, in particular along the optical disk of M31 and M33. 
We expect that most of these sources are actually candidates asymptotic giant branch 
(AGB) stars, planetary nebulae (PNe), young stellar objects (YSO) or red stars. These 
are interesting objects to study in their own right. Further more, finding quasars in the 
very central regions of M31 and M33 will be extremely valuable for the PM measurements 
and ISM spectroscopic studies. As such we have included all those candidates in our input 
catalogs for LAMOST observations. The cyan and blue circles in Fig.[T] represent the 
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LAMOST fields observed during the 2010 and 2011 observational seasons; the sky area 
covered during each year was ~ 100 deg^ and ~ 135 deg^, respectively. In 2010, only 
optically selected low-redshift quasars candidates based on the SDSS colors were selected. 

3. Observations and Data Reduction 

M31 has been targeted by the LAMOST since the commissioning phase initiated in 
2009. Paper I presented some of the early results about background quasars discovered 
by LAMOST. Here we present results based on data collected in the 2010 and 2011 
observational seasons. 

The exposure time for the plates vary from 600 to 1800s, repeated 2 or 3 times to 
allow rejection of cosmic rays and to increase the signal-to-noise ratio (S/N). In total, 731 
and 5641 unique quasar candidates were observed in 2010 and 2011, respectively. Some 
sources were targeted repeatedly. This is partly because the field of view of LAMOST 
is circular, which forces adjacent fields to overlap to allow coverage of a large contiguous 
sky area. In 2010, only SDSS optically-selected low-redshift quasars candidates were 
selected and observed. In October and November of 2010, LAMOST CCDs suffered severe 
electromagnetic interference and insufficient dewar cooling, and so few useful data were 
collected. See Fig.[I]for the fields targeted by the LAMOST in both 2010 and 2011. 

The LAMOST has 16 low resolution spectrographs, each accommodating 250 fibers. 
From the fall of 2011, LAMOST impleted slit masks with a width of 2/3 the 3.3arcsec 
fiber diameter, yielding a spectral resolving power ~ 1800. Before that a mixture of full 
(i.e. no slit masks) and 1/2 slit modes were used, but mostly in full slit mode, which yields 
a spectral resolution ~ 1000. The light that enter each spectrograph is split into two 
channels, the blue and red channels, for which the wavelength coverages are 3700-5900 A 
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Fig. 1. — Fields observed with the LAMOST in the vicinity of M 31 and M 33 during the 2010 
and 2011 observational seasons. The positions of M31 and M33 are labeled with characters, 
M31 and M33 are centered in the white areas. The grey area represents the SDSS footprint 
in this area. The magenta plus symbols delineate the centers of the individual XSTPS fields 
covered the M31/M33 area, with 2x2 sq. degrees for each field. The cyan and blue circles 
represent the LAMOST fields observed in 2010 and 2011, respectively. The black circles 
represent the three SDSS spectroscopic fields targeting quasar candidates, the two plates in 
the northeast halo have the same central positions. 
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and 5700-9000 A, respectively. The spectra from each spectrograph are imaged onto two 
4096x4096 CCD cameras (See Cui et al. 2012 for further details). The spectra were reduced 
using the LAMOST 2D pipeline, including bias subtraction, spectra tracing and extraction, 
wavelength calibration, flat-fielding and sky subtraction (see Luo, Zhang & Zhao 2004; 
Luo et al. 2012). Considering the uncertainties in the commissioning phase and the pilot 
surveys, and the relatively low S/N of most quasar candidate spectra, for the current work 
of identifying quasars, the non flux-calibrated spectra were used. Hence the photometric 
magnitudes serve as much better measurements of the brightness of our targets. 

4. Results and Discussion 

This current work is based on observations collected during the LAMOST commissioning 
phase and the pilot surveys when the LAMOST performance was still being characterized 
and optimized. Thus it is unavoidable that the data presented here suffer from some defects 
in one way or another, such as low throughput, poor sky subtraction and so on. The 
fiber positioning, although improved significantly in spring 2011, still has room for further 
improvement. Fortunately, given their characteristic broad emission line, quasars are easily 
identified. We visually examined the ID spectra of each individual quasar candidate. We 
required that at least two emission lines were securely identified. The results led to the 
discovery of 17 and 509 new quasars in the vicinity of M31 and M33 in the fields observed 
in 2010 and 2011, respectively. If we adopt a luminosity cut following the SDSS Quasar 
Survey (Schneider et al. 2010, see also the references therein), with absolute magnitude 
Mi = —22.0 in a cosmology H q = 70 km s^^ Mpc^^, Qm = 0.3, and ^a = 0.7, then the 
number of newly-discovered quasars are 15 and 509 in 2010 and 2011 datasets, respectively. 

All the new quasars have reasonable S/N and at least two identified emission lines, 
thereby allowing reliable redshift estimates (see Fig.|2]). In Tabled! we present the catalog 
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Fig. 2. — Example of LAMOST spectra of newly-discovered quasars, with identified lines 
labeled. The relative flux is in units of counts per pixel, the horizontal dotted line labeled 
the zero flux. The spectra were binned by a factor of 4, and remaining cosmic rays after 
pipeline processing were removed manually for clarity. 
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Table 1. Catalog of New Quasars in the Vicinity of M31 and M33 Discovered by 
LAMOST in Fields Observed in 2010 and 2011. 



Object 


R.A. 


Dec. (J2000) 


Rcdshift 


u 


9 


r 


i 


Z 


A(i) 


Selection 


J003255. 71+394619.2 


8.232129 


39.772018 


1.134 


_ 


18.53 


18.24 


18.29 


_ 


0.12 


w,-,- 


J003432. 52+391836.1 


8.635497 


39.310022 


0.138 


18.55 


18.53 


18.38 


17.87 


18.03 


0.10 


-,S,M 


J003459. 61+420655.2 


8.748376 


42.115341 


2.472 


- 


18.86 


18.60 


18.71 


- 


0.12 


w,-,- 


J003506. 69+404003.4 


8.777890 


40.667629 


1.844 


20.27 


20.05 


19.91 


19.56 


19.43 


0.14 


-,s,- 


J003514.29+401414.6 


8.809581 


40.237413 


0.279 


- 


19.50 


18.74 


18.57 


- 


0.12 


w,-,- 


J003522.68+393353.9 


8.844538 


39.564987 


0.670 


- 


19.32 


19.48 


19.53 


- 


0.13 


w,-,- 


J003524.49+394619.1 


8.852052 


39.771992 


2.340 


- 


19.25 


19.03 


18.77 


- 


0.12 


w,-,- 


J003535. 01+404351. 2 


8.895891 


40.730891 


1.774 


20.33 


20.21 


20.25 


19.90 


20.12 


0.13 


-,s,- 


J003546. 45+413358.6 


8.943565 


41.566297 


2.580 


- 


19.21 


18.93 


19.02 


- 


0.14 


w,-,- 


J003619. 34+411746.5 


9.080616 


41.296262 


1.651 


20.38 


19.90 


19.56 


19.21 


19.09 


0.14 


-,s,- 



Note. — W represents candidates selected by the IR criteria based on the WISE data, S represents those optically selected with 
the SDSS criteria for low-redshift quasars, M represents those optically selected from the KPNO 4m Local Group Galaxy survey 
data (Massey et al. 2006). Only a portion of the Table is shown here for illustration. The whole Table contains information about 
526 new quasars, and is available online in the electronic version. 
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of these newly discovered quasars. The properties of these quasars, including target names 
(in format of Jhhmmss.ss+ddmmss.s), equitorial coordinates, redshifts, the observed SDSS 
u, g, r, i, z magnitudes without extinction corrections, z-band extinction from Schlegel 
et al. (1998), and the selection criteria, are listed. Only a portion of the Table is shown 
here. The whole Table, containing information about all 526 new quasars, is available in 
the online electronic version. In the last column of the Tabledl W indicates candidates 
selected with the IR criteria based on the WISE data, S indicates those optically selected 
using the SDSS criteria for low-redshift quasars, and M indicates those optically selected 
candidates based on the data from the KPNO survey of Local Group Galaxies. Columns 
5-9 give the SDSS u, g, r, i, z magnitudes. For some targets only the g^ r, i magnitudes 
are given. Those are targets selected by cross-correlating the WISE-IR candidates with 
the XSTPS optical catalog. Fig. [2] shows spectra of 7 new quasars spanning a wide range 
of redshifts. The spectra have not been flux-calibrated. At wavelengths shorter than 4000 
A, where the instrument throughputs are low, as well as between the dichroic cross-over 
wavelength range 5700-6000 A, large spike- and trough-like artifacts often appear in the 
flux-calibrated spectra. This is particularly true for spectra with relatively low S/N, as is 
the case for many of our targets of interest here. Because of this, we have opted to not to 
use the flux-calibrated spectra. Since the LAMOST spectra were oversampled, we binned 
the spectra by a factor of 4 to improve the S/N. Cosmic rays which remained after pipeline 
processing were removed manually for clarity. The catalog and the individual spectra are 



available at the LAMOST public website http://www.lamost.org 



In addition to the 526 conflrmed quasars, another 20 and 193 'probable' quasars 
have been identifled in the current LAMOST 2010 and 2011 datasets. They either have 
marginal S/N or have only one emission line detected and so are at uncertain redshifts. 
Further observations are needed to conflrm their identiflcations. We take this 'probable' 
quasar catalog as a supplemental list, this supplemental list and the related spectra are also 
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available at the aforementioned website, which will be updated in real time as new data 
are obtained or the pipeline is improved. Finally, 11 and 23 previously known quasars are 
re-observed in the LAMOST 2010 and 2011 datasets, respectively. 

The sample has a magnitude limit of z = 20.0. Amongst the 526 newly discovered 
quasars, the brightest, J010415. 77+402243. 9, has an i magnitude 15.45 with redshift of 
0.193, and was previously unknown (its spectrum is shown in Fig.|2]). There are 5/20/75 new 
quasars brighter than 17.0/17.5/18.0 magnitude in z-band, respectively. For comparison, 
amongst the quasars discovered in the three SDSS spectroscopic plates (Adelman-McCarthy 
et al. 2006, 2007) and in the LAMOST 2009 commissioning observations, the corresponding 
numbers are 1/1/12. Within a lOdeg. radius of M31 and of M33, there are 155 previously 
known quasars with redshift estimates reported in the NED archive (see also Fig.[6]). 
Many of them have multiple photometric measurements taken in a variety of filters, 
and 151 of them fall in the area covered by the SDSS or XSTPS surveys. Using a 
3" cross-correlation radius, we searched their counterparts in the SDSS and XSTPS 
sources with i magnitudes brighter than 21.0 and 20.0, respectively. We find 116 of 
them have i band counterparts, amongst them there are 20/41/52 quasars brighter than 
17.0/17.5/18.0 magnitudes, respectively. These bright quasars will serve as good candidates 
for absorption-line spectroscopic studies of the ISM of M31, M33 and of their related 
substructures. Fig.[3] shows the magnitude distributions of the 532 quasars detected in the 
2011 datasets, including 509 newly discovered quasars and 23 previously known sources 
which were re-observed by LAMOST. The magnitude has a bin size of 0.1, the y-axis is in 
quasar number density. In the magnitude distributions, as well as the following redshift 
distributions, only quasars detected in ~ 135 sq. deg in 2011 are included, the quasars 
detected in 2010 are not considered. The magnitude distributions of the SDSS DR7 quasars 
(Schneider et al. 2010) are also shown for comparison. We can see that a large number of 
quasars with i magnitude fainter than 18.0 are left out by LAMOST due to its relatively 
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Fig. 3. — i-band magnitude distributions of the 532 quasars detected by LAMOST in 2011 
(grey), including 509 newly discovered and 23 previously known quasars detected by LAM- 
OST. The i-band magnitude distributions of the SDSS DR7 quasars are shown here for 
comparison (black), the steep gradients at i ps 19.1 and i ^ 20.2 are due to the magni- 
tude limits of the SDSS selection algorithms for the low- and high-redshift quasars. The 
magnitude bin size is set to 0.1. The y-axis is displayed in quasar number density. 
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low performance during the commissioning and pilot survey phases. 

Fig.HJshows the redshift distributions of the 532 quasars detected by LAMOST in 2011. 
The redshift has a bin size of 0.1, the y-axis is in quasar number density, as in Fig.O The 
lowest redshift is found for J003842.80+384551.9, at ^ = 0.093, whilst the highest redshift 
object is J011118. 08+325833. 8, a.t z = 3.204. For the optical selection criteria adopted here, 
low-redshift quasars selected as outliers of stellar loci on the SDSS color-color diagrams 
are mainly z < 2.2 quasars (Fan 1999; Richards et al. 2001), whereas the optical/IR 
algorithm is capable of finding quasars of redshifts up to z < 3.5 (Wu et al. 2012). We 
notice that there is no significant difference in redshift distribution between the optical and 
optical/IR selected quasars, as illustrated by the blue and red lines in Fig.HJ respectively. 
For comparison, we also plot the redshift distributions of SDSS DR7 quasars. Except for 
the relatively lower number-density in redshift space, we can also clearly see that there are 
two troughs in the redshift distributions of LAMOST detected quasars (see Fig.H]), one 
near redshift ~ 0.9 and another near ~ 2.1. These are entirely due to the selection effects. 
For redshift z < 0.9, the Mg II A2800, B.6, H7, H/3 and Ha, as well as the forbidden lines 
[O III] A4959 and A5007 all fall within the LAMOST wavelength coverage. For redshifts 
l.l < z < 2.0, the C IV A1549, C ill] A1908, and Mg ll A2800 lines are redshifted into the 
LAMOST wavelength range. For redshifts z > 2.1, Lya and N v. Si IV and C iv resonant 
lines as well as the C ill] A1908 line become visible. For quasars near the z ~ 0.9 trough, 
only the Mg II A2800 line is easily detectable, other emission lines such as B.6, H7 and 
H/3 are all redshifted to wavelengths longer than 7000 A, where the spectra are seriously 
contaminated by telluric emission lines. The difficulty remains until the redshift reaches 
z ~ 1.1 when the C III] A1908 becomes visible. At redshifts ~ 2.1, although there are 
three emission lines, the C iv A1549, C III] A1908, and Mg ll A2800 lines fall in the optical, 
the C III] A1908 line is unfortunately redshifted to a wavelength of 5900 A, the cross-over 
wavelength of the dichroic where the instrument has a very low throughput, and the Mg ll 
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Fig. 4. — Redshift distributions of the 532 LAMOST detected quasars in 2011 (thick grey). 
The blue and red hues represent redshift distributions of SDSS low-redshift and optical/IR 
selected quasars, respectively. The redshift distributions of the SDSS DR7 quasars are also 
shown here for comparison (black), the low number density at 2; ~ 2.7 is caused by the 
degeneracy of the SDSS color of stars and colors of quasars near this redshift. The redshift 
bin size is set to 0.1. The y-axis is displayed in quasar number density. 
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A2800 is redshifted to a wavelength ~ 8700 A where sky emission hnes dominate the region. 
For the 213 'probable' quasar candidates detected in the LAMOST 2010 and 2011 datasets, 
either only one emission line had been detected or the spectrum was of marginal S/N. Since 
we were unable to discriminate the C IV A1549 and Mg ll A2800 lines, no redshifts were 
provided. Fig.[5] shows the distributions of quasars in redshift-magnitude space, we can see 
that not only the selection effects in redshift space, but also the low number density of 
quasars for z-band fainter than 18.0, indicate that the performance of LAMOST and 2D 
pipeline need to be carefully characterized and optimized. 

Fig.ini plots the spatial distribution in the ^-r] plane of all known background quasars 
found in the vicinity of M 31 and M 33, including the 526 newly identified with the LAMOST 
which are reported in this paper, 14 quasars discovered by the LAMOST based on earlier 
commissioning observations (Paper I), 75 SDSS quasars (Adelman-McCarthy et al. 2006, 
2007) and 155 previously known quasars with redshifts listed in the NED archive. Here ^ 
and r] are respectively Right Ascension and Declination offsets relative to the optical center 
of M31 (Huchra, Brodie & Kent 1991). In Fig.Ej the magenta stars denote the central 
positions of M31 and M33, respectively. The magenta ellipse represents the optical disk of 
M31 with an optical radius R25 = 95'3 (de Vaucouleurs et al. 1991), an inclination angle 
i = 77deg. and a position angle PA = 35deg. (Walterbos & Kennicutt 1987). The quasars 
are distributed around M31 and M33, along the Giant Stellar Stream (McConnachie et 
al. 2009), as well as the extended halo of M31. Amongst them, there are 93 quasars 
within 2.5° of M31 (about 34kpc assuming a distance of 785 kpc to M31; McConnachie 
et al. 2005), 73 of them are newly discovered, 7 are from Paper I and 13 are previously 
known quasars listed in the NED. Several quasars fall very close to the optical disk of 
M31. Tens of quasars behind the Giant Stellar Stream have been found, and the number 
of known quasars behind M31's extended halo and its associated substructures increases 
by a substantial amount. The much-enlarged number of known quasars in the vicinity of 



-21 - 




Fig. 5. — Redshift versus z-band magnitude distributions of the 532 LAMOST detected 
quasars in 2011 (blue points). The contours represent the redshift versus z-band magnitude 
distributions of the SDSS DR7 quasars. The steep gradients at i ~ 19.1 and i ~ 20.2 are 
due to the magnitude hmits of the SDSS selection algorithms for the low- and high-redshift 
quasars; the low number density at 2; ~ 2.7 is caused by the degeneracy of the SDSS color 
of stars and colors of quasars near this redshift. 
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M31 and M33 should provide an invaluable sample for future PM and ISM/IGM studies of 
M31, M33 and the Local Group. 

No quasars have been identified near the central regions of M31 and M33 due to 
crowding in the field. Outside those regions, the spatial distributions of identified quasars 
is not uniform in the 11 fields observed in 2011, one covering M33 yielded 147 quasars, 
giving the highest quasar number density of ~7.5/sq. deg. On the other hand, a plate 
near the bottom-left corner of Fig. [6] yielded only 34 quasars. The quasar number density 
identified by the LAMOST is lower than that of SDSS for low-redshift quasars with 
point-like morphology, approximately 10 deg"^ (Schneider et al. 2010). The low yields 
were mainly due to insufficient S/N, either because of poor weather conditions or lower 
than expected performance of the LAMOST. Since the current sample of quasars were 
identified using the commissioning and pilot survey data, the quasar selection algorithm 
and telescope performance was still undergoing development and testing, hence the sample 
is not homogeneous and is not intended for statistical analysis. 

This work is partially supported by the Young Researcher Grant of National 
Astronomical Observatories, Chinese Academy of Sciences, and China Postdoctoral Science 
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Fig. 6. — Spatial distributions of background quasars in the vicinity of M31 and M33. 
Orange, cyan and blue filled circles represent quasars identified in the LAMOST 2009, 2010 
and 2011 datasets, respectively. Crosses and open squares represent SDSS quasars and 
previously known quasars with redshifts archived in the NED, respectively. The magenta 
stars mark the central positions of M31 and M33, while the magenta ellipse represents the 
optical disk of M31 of radius R25 = 95'3. 
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